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Abstract N-body simulations predict that dark matter halos with different mass scales 
are described by a universal model, the Navarro-Frenk- White (NFW) density profiles. As 
a consequence of baryonic cooling effects, the halos will become more concentrated, and 
similar to an isothermal sphere over large range in radii (~ 300 /i~ 1 kpc). The singular 
isothermal sphere model however has to be truncated artificially at large radii since it 
extends to infinity. We model a massive galaxy halo as a combination of an isothermal 
sphere and an NFW density profile. We give an approximation for the mass concentration 
at different baryon fractions and present exact expressions for the weak lensing shear and 
flexion for such a halo. We compare the lensing properties with a Singular Isothermal 
Sphere and NFW profiles. We find that the combined profile can generate higher order 
lensing signals at small radii and is more efficient in generating strong lensing events. 
In order to distinguish such a halo profile from the SIS or NFW profiles, one needs to 
combine strong and weak lensing constraints on small and large radii. 
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1 INTRODUCTION 

The Cold Dark Matter with the cosmological constant model (ACDM) provides a successful descrip- 
tion of many properties of observations of the universe. N-body simulations of ACDM models predict 
dark matter halos with a universal density profile (e.g. Navarro et al., 1997). The Navarro-Frenk- White 
(NFW) profile appears to be a good approximation for dark halo profiles over a wide range of masses. 
On the other hand, the NFW halo density profile can also be generalized with an arbitrary power law 
central cusp, and outer regions that fall off as ?'~ 3 (Jing & Suto, 2000). It also has been found that the 
slope of the inner regions steepens for smaller mass haloes. More importantly, baryonic cooling will sig- 
nificantly steepen the density profiles, close to the isothermal slopes observed (Koopmans et al., 2009). 
The baryon effect is more significant in the galaxy halo since it contains more baryons. A composite 
model with an NFW dark matter halo and a de Vaucouleurs stellar component is suggested for mas- 
sive galaxies by Gavazzi et al. (2007). The total density profile is close to isothermal form over large 
range in radius (~ 300 /i _1 kpc). Therefore, we model the halo total mass profile as an Isothermal-NFW 
(INFW) profile, which is the combination of an NFW dark halo plus a stellar component at inner radii, 
i.e. p oc r~~ 2 for small radius. 

Gravitational lensing provides a direct way to study the mass distribution of large scale structures in 
the universe as well as galaxy and cluster halos. It probes the mass distribution independent of the nature 
of matter or its dynamical state (e.g. Bartelmann & Schneider, 2001; Treu, 2010). Lensing is widely 
used for the cluster mass reconstruction (e.g. Bradac et al., 2006), and galaxy halo measurement (e.g. 
Cacciato et al., 2009). Weak lensing is the physical phenomenon causing the weak image distortion of 
background galaxies. By comparing the image distortions with non-lensed image shapes, one can infer 
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the mass distribution of the foreground lens. In weak lensing, most studies consider the shear effect, 
which transfers a round source into an elliptical one. Higher order effects, flexion are gradually coming 
within reach. Flexion can be introduced as derivatives of either the surface mass density or the shear. 
They respond to smaller-scale variations in the projected mass distribution than the shear (Bacon et al., 
2006). The convergence gradient, called the first flexion T, introduces a centroid shift in the lensed 
image and is a spin-1 symmetry quantity, while the second flexion Q is the gradient of shear and is spin- 
3. Flexion provides a measure of small scale variations of mass distribution as well as the halo ellipticity 
(Er & Schneider, 201 1; Er et al., 2012). 

The lensing properties of different halo profiles have been widely studied, e.g. the NFW profile 
(Bartelmann, 1996) and the Einasto profile (Re tana-Montenegro et al., 2012). Wyithe et al. (2001) and 
Keeton & Madau (2001) have studied a generalized NFW type profile for lensing. Therefore, it is in- 
teresting to use the INFW profile as a galaxy halo, and it is natural to extend its applications to the 
gravitational lensing characteristics of dark matter halos. For first time, we apply analytical and numer- 
ical methods to the gravitational lensing study of INFW halo profiles. In Sect.2, we present the basic 
halo properties of the INFW profile. In Sect.3, the analytical formula of an INFW lens halo is given. 
We compare the INFW profile with other models in Sect.4 and give a summary at the end. The cosmol- 
ogy that we adopt in this paper is a ACDM model with parameters based on the results of the Wilkinson 
Microwave Anisotropy Probe seven year data (Komatsu et al., 201 1): Oa = 0.734, O m = 0.266, Hubble 
constant H = lOO/i kms -1 Mpc -1 and h = 0.71. 



2 INFW HALO PROPERTIES 



In analogy to NFW model, the density profile of INFW is given by 



p(r) 



r 2 (r s + r) 



(1) 



where p c = 2>H (z) 2 / (8irG) is the critical density of the universe, H(z) is the Hubble parameter, and G 
is Newton's constant. The dimensionless characteristic density is given by 
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(see, e.g., Wyithe et al., 2001). We will use the same definition for the concentration cj = r2oo/ r s, 
where r s is the scale radius. The virial radius r 2 oo is defined as the radius inside which the mass density 
of the halo is equal to 200p c (Navarro et al., 1997). The mass of a halo contained within a radius of r2oo 
is thus 
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There is no specific study or simulations for the relationship between mass and concentration for 
generalized-NFW profiles. We assume that initially dark matter and baryons follow the same NFW 
profile. Due to the cooling effect, baryons collapse toward the center of the dark matter halo, and steepen 
the inner density profile. We assume the collapsed baryons make a fraction of the total mass. The 
outer density of the INFW profile will become lower by a factor of (1 — /b): /?infw ( r 200 ) = (1 — 
./6)PNFw( r 20o)- We take ft as the universal baryon fraction, although a lower number does not change 
the scaling significantly. A relation between c and c/ can be obtained from 
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This relation can be solved numerically. It can be also approximated by 
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Fig. 1 The approximate relationship between the concentration c for the NFW profile and a for the 
INFW profile. The points are the numerical results from solving Eq.(4) for different baryon fractions: 
fb = 0.16 (circles), ft = 0.08 (squares) and fb = 0.06 (triangles). The lines are the approximate 
relationship using Eq.(5). 




In Fig. 1, one can see that our approximation mainly agrees with the numerical results. A smaller baryon 
fraction will lead to a lower concentration of the INFW halo. When the concentration ci becomes to 0, 
r s —> oo, and the INFW profile reduces to an SIS. Thus in general the INFW profile is more concentrated 
than SIS profile at small radial. We will see in next section that the INFW profile can produce higher 
lensing signals and is more efficient in generating strong lensing than other profiles. The small variation 
of baryon fraction does not strongly affect the matter density profile (right panel of Fig. 2). With higher 
baryon fraction, the density at inner radius is larger. In the rest of the paper, we will use ft = 0.16 and 
Eq.(5)to estimate the concentration of the INFW halo. 

In left panel of Fig. 2 we show p(r) for three different halo profiles using same halo mass M200 One 
can see that the INFW profile has the same slope as SIS at small radii (< 30 hr 1 kpc) and approaches 
to NFW at large radii. 

3 LENSING PROPERTIES OF INFW HALO 
3.1 Basic lensing formula 

The fundamentals of gravitational lensing can be found in Bartelmann & Schneider (2001). For its ele- 
gance and brevity, we shall use the complex notation. The thin-lens approximation is adopted, implying 
that the lensing mass distribution can be projected onto the lens plane perpendicular to the line-of-sight. 
We introduce angular coordinates 9 on the lensing plane with respect to the line-of-sight. The lensing 
convergence, that is the dimensionless projected surface-mass density, can be written as 

K(0) = £(0)/£ cr , where E cr = ° D ' (6) 

47T(_r i^d^ds 

is the critical surface mass density depending on the angular-diameter distances D s , Z?d and Z?d s from 
the observer to the source, the observer to the lens, and the lens to the source, respectively. £(0) is 
the projected surface-mass density of the lens. All lensing quantities can be derived from the effective 
lensing potential ip, 

il>{6) = \\ & 2 Q'k{0') ln|0 - 0'\ . (7) 
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The lens equation is given by 



P = e-a(6), (8) 
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Fig. 2 Left panel: halo mass density p(f) for three different profiles: NFW (solid line), INFW (dashed 
line), SIS (dot-dashed line). The same mass (A/200 = 10 12 hT 1 M@) is used for different profiles (also 
for right panel). The concentration is c = 6.95 (cr = 3.11, f b = 0.16) for NFW (INFW) halo. Right 
panel: p(r)iNFw/p(r)sis w ^ tn different /;,: 0.16 (solid line), 0.08 (dashed line), 0.06 (dot-dashed line). 



where j3 is the source position and a is the deflection angle 

a = V c ^, (9) 
where the complex differential operators is defined as 

d d , d d 

To the lowest order, image distortions caused by gravitational lensing are described by the complex 
shear and convergence (which equals to the dimensionless surface mass density) 

7 = \ (%4> - + idxd^ = ^V^; k = \{d\^ + 9 2 V) = ^V C V^, (11) 

where the subscripts i denote partial derivatives with respect to 9i. The magnification for a point source 
is thus given by 

(1 - K y - I7I 2 

The shear transforms a hypothetical round source into an elliptical image. The T and Q flexions can be 
introduced as the complex derivatives 

T = V c k; Q = V c7 . (13) 

The flexions are thus combinations of third-order derivatives of the effective lensing potential ip. We 
shall denote their real and imaginary parts by (J 7 , Q)\ and (J 7 , 0)2, respectively. In terms of the lensing 
potential, we have 

T = T x + LF 2 = 1 (Oft/j + dtd^) + i (dld 2 i> + d^j) (14) 



and 



G = 01+102 = 1 {dli> - Zdrdli>) + i (3dfd 2 ^ - dlif) . (15) 
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3.2 Lensing of INFW halo 



We derive the analytical expression for the lensing properties of INFW halo. The surface mass density 
of a spherically symmetric lens is obtained by integrating along the line of sight of the three-dimensional 
density profile, 



s(0 



pWe + z 2 )dz, 



(16) 



where £ is the distance from the center of the lens in the projected lens plane £ = 9D S . It implies the 
following form for the dimensionless surface mass density 



k(x) = 2k s (J^- - f(x)J 



(17) 



where x = 9/8 s [0 S = r s /D^j, and f(x) is given by 
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In the spherical case, the deflection angle is given by 
2 
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The analytical form of the shear can be calculated from j(8) = [k(9) — k(9)} cxp[2i0], where <f> is the 
polar angle. R{6) is the mean surface mass density within a circle of radius of 9 from the lens center (see 
e.g. Bartelmann & Schneider, 2001). The expression for shear due to the INFW is 



7(.t) = 2k s 
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where f{x) is defined in Eq. (18). The analytical form of two flexions can be also given 
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The elliptical INFW lensing properties can be calculated numerically (Keeton, 2001). 

Furthermore, as pointed out by Schneider & Seitz (1995); Schneider & Er (2008), due to the mass- 
sheet degeneracy, the directly measurable properties are the reduced shear and reduced flexion 
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The weak lensing properties of the INFW profile also show approximated behavior as a combination 
of two power-law profiles. At small radii, the asymptotic behavior can be approximated by a SIS, i.e. 



k, 7 oc u , and J 7 , G oc 9 2 . At large radii, it behaves like the power-law p oc 
signal rapidly fade out, k, 7 oc 9~ 2 , and J 7 , G oc 8~ 3 . 



Thus the lensing 
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4 PROFILES COMPARISONS 

We compare the weak lensing properties for INFW, NFW, and SIS profiles. We use an approximation 
relation to calculate the mass concentration of the NFW profile (Neto et al., 2007) 



and use Eq.5 to obtain cj. The velocity dispersion of the SIS profile a v is calculated through a\ = 
-j^GH(z) M 2 oo (Mo et al., 1998). The lensing properties of the SIS or NFW profiles can be found in 

e.g. Wright & Brainerd (2000); Schneider et al. (2006). We use lens halo mass M 20 o = 10 12 hr 1 Mq, 
which is a galactic sized halo. The lens is placed at redshift Zd = 0.2 and the sources are at redshift z s = 
1.0, which are accessible median redshifts for galaxy survey, e.g. SDSS or LSST The concentration 
parameter for the NFW (INFW) profile in our test is c = 6.95 (c/ = 3.11). The Einstein radius of the 
SIS profile is 6>e ~ 0.3 arcsec. 

Fig. 3 shows the predicted convergence, reduced shear, first and second reduced flexions as a func- 
tion of the angular separation from the lens center. The mass profile of the mock galaxy halo is assumed 
to be SIS (dotted line), NFW (dot-dashed line) and INFW (solid line) model. One can see that the overall 
behaviors of the three profiles are comparable. The asymptotic lensing behavior of the INFW profile are 
proportional to that of the SIS profile at small radii and approach NFW profile at large radii. The signal 
magnitudes of all lensing properties for the INFW halo are stronger than the other two at small radii but 
drop faster and eventually below that of the other two profiles. The differences between the magnitudes 
of the lensing signal are stronger at small radii than that at large radii. In particular, the shear and second 
flexion show a great dissimilarity. At large radii, the difference between three profiles is not significant. 

In particular, the weak lensing properties of NFW halos is significantly different from the other 
two profiles at small radius (< 10 arcsec). The lensing signal of NFW is shallower than the others. On 
the other hand, the signal magnitude of INFW halo is about 2 times higher than that of SIS halo. The 
lensing signals of INFW drop faster with increasing radius, as one can see clearly from the first flexion 
(the convergence k is not an observable quantity). Moreover, the angular separation at which the INFW 
halo first flexion is exceeded by other profiles is around 10 arcsec, and it is larger for shear and second 
flexion. In principle one can study the weak lensing signal, i.e. the shear and flexion to constrain the 
halo density profile. However, the weak lensing signal at large radii is small and hard to detect. On the 
other hand, it is also difficult to measure weak lensing signal when the background image is close to the 
lens galaxy. One can perform stacking method for galaxy-galaxy lensing studies. A large volume survey 
is necessary. 

Moreover, the significant difference lensing properties at small radius will cause different strong 
lensing signal. In order to simply see the strong lensing properties, we compare the magnification cross 
sections for the three profiles. The cross section for a given magnification threshold is defined as 



In Fig. 4, the magnification cross section is shown in the left (right) panel for halo with mass 
10 13 ft -1 Mq (different halo mass). We can see from both panels that the INFW profile can gener- 
ate larger cross sections than the other two profiles, due to the high mass concentration of the INFW 
profile (top left panel in Fig. 3). The cross section of NFW halo increases faster with mass than other 
profiles, but decrease faster with [i. The curves of INFW and SIS profiles again have similar shapes, 
but the cross section of INFW halo is about two times larger than that of SIS halo for halo mass of 
10 13 h~ Y Mq. In additional tests, we also study the cross section of strong lensing multiple images. The 
probability generated by INFW halo can be several times higher than NFW halo, and will be easy to 
distinguish from each other. On the other hand, the INFW model generates about 3 times higher multi- 
ple image cross section than SIS model with halo mass ~ 10 12 h^ 1 Mq and approaches to that of SIS 
model for massive halo (> 10 15 hr 1 Mq). The concentration cj becomes small for massive halo, thus 
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Fig. 3 Convergence, reduced shear and reduced flexions of three different profiles: INFW (solid line), 
SIS (dashed line), NFW (dot-dashed line). The mass of the lens halo is M200 = 10 12 hT 1 Mq. The lens 
and source redshifts are assumed to be Zd = 0.2 and z 3 = 1.0 respectively. 



the INFW nearly reduces to SIS profile. The multiple image separation generated by INFW lens can 
reach 4 arcsec for a halo mass of 10 13 ft -1 Mq, which is about 40 percent larger than that generated 
by SIS lens. Therefore, the galaxy-galaxy strong lensing statistics can be a potential tool to distinguish 
INFW and SIS profile. 

5 SUMMARY 

We have studied the lensing properties of the INFW mass profile. The INFW profile is motivated by the 
combination of Cold Dark Matter simulations and a stellar component in the inner region of the dark 
matter halo, together with some evidence from observations (Gavazzi et al., 2007). The inner profile of 
INFW is isothermal, i.e. p oc r~ 2 and the outer profile is NFW-like p oc r~ 3 . An approximate mass 
concentration due to different baryon fractions is given for the INFW profile, as a direct consequence of 
baryon collapse toward the center of halo. 

The analytical expressions for deflection angle, convergence, shear and flexions of an INFW halo 
lens are given. We have compared the lensing properties of INFW profile with NFW and SIS halo pro- 
files. We find that the INFW profile is more efficient than the others in generating lensing magnification, 
and the weak lensing signals of INFW halo is stronger at small radii than that of other profiles for the 
same halo mass. Strong lensing statistics can be used to constrain the lens profile, e.g. the image separa- 
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Fig. 4 The cross section of lensing magnification cr(fi) (in unit of arcsec 2 ) for three halo profiles: 
INFW (solid line), SIS (dashed line), NFW (dot-dashed line). The left panel shows the cross section as 
a function of lensing magnification for a lens halo with mass M200 = 10 13 h^ 1 Mq. The right panel 
shows the cross section of given magnification (/1 = 2) for different lens halo mass. Same redshift 
condition of Fig.3 (zd = 0.2, z s — 1.0) is used. 



tion. However, the image separation statistics is only sensitive to the inner profile of the lens halo. There 
is a degeneracy between the massive SIS lens and high concentrated INFW lens. Therefore, using weak 
lensing to study the large radial profile is essentially necessary. 
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